SUNY-NTG-97-50 LA-UR-97-4446 



Instanton-induced charm contribution 
to polarized deep- inelastic scattering 

00 ' 

Q\ ', Andree Blotz W and Edward Shuryak ^ 

(1 ) Theoretical Division, T5 MS B283, LANL, Los Alamos, NM 87545, USA 

X> , 

r ^ ! email:blotz@t5.lanl.qov 

. 

(2) Department of Physics /Astronomy, State University of New York at Stony Brook, 
^ ! NY 11794-3800, USA 

> : 

■^j- , email:shuryak@dau. physics. sunysb.edu 

g ; (February 7, 2008) 



ON 

a: 



Abstract 



^ e— I Recent data on B decays involving n' may be explained if the singlet 

axial current of charmed quarks has a large matrix element to r]', and 

x. 

instantons were shown to be able to generate this effect. We study the 
magnitude of charm contributions to nucleon polarized structure func- 
tions generated in a similar way. Comparing the charm contribution, 
which is related to a dim(6) gluonic operator, to that of light quarks, 
which are related by the anomaly equation to GG, we found that Ac/AS 
= -(0.2-0.08). Future experiments like COMPASS at CERN identifying 
charm production in DIS can measure this "intrinsic polarized charm" 
component of the nucleon. 
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1. It was believed that non-perturbative phenomena in QCD can generate light 
quark (and gluon) "sea" in nucleon structure functions, while the admixture of heavy 
flavors such as charm is (i) very small due to their mass, and that (ii) it should 
be produced by perturbative effects alone. Although the first part of the argument 
(suppression by the power of charm mass squared) is correct, it can be counterbalanced 
by sufficiently strong non perturbative gluon fields. Small-size instantons present in the 
QCD vacuum have fields G ~ lGeV 2 , which is not much smaller than m? c ~ 2GeV 2 , 
and so in principle one may expect a sizable "intrinsic charm" of at least some hadrons 
(as anticipated by Brodsky et al. for a long time [T|]). 

Probably the first indication of that came from recent CLEO data on B — > rj'K 
and inclusive B — > 77' + . . . decays @ , which can be explained H if r/' has such a charm 
component. Quantitative estimates of this effect - based on an instanton liquid model 
for the non perturbative vacuum - were made in Q, which concluded that indeed 
current information about instantons lead to surprisingly large charm content of the 
7/, comparable to what is needed to explain the CLEO data. 

It was further speculated in |4| that other processes related to axial current of 
charm quarks may also show similar effects. One of them is the charm contribution 
to the polarized deep inelastic scattering (DIS), which we discuss in this work. Crude 
estimates given in 0] suggested that it may be even comparable to that of all light 
quarks together (which is, we remind, already strongly suppressed, as compared to the 
naive unit value for the "spin content" of the nucleon). 

In the present paper we address the question of an intrinsic polarized charm com- 
ponent directly, by evaluating relevant three-point correlation functions in the vacuum. 
We use the same instanton liquid models as in ref. 0. The results however are about 
an order of magnitude smaller than the estimates presented there. Nevertheless, they 
are still very different from the perturbative predictions, which imply that the charm 
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sea is mostly unpolarized. 

2. Let us first recall few standard formulae for polarized DIS on a nucleon. [] The 
first moment of the polarized structure function 

T?(Q 2 )= ( l dx 9l (x,Q 2 ) (1) 
Jo 

can be related to the matrix element of the axial current by the standard OPE treat- 
ment. The answer (including charm component Ac) can be written as ffl 



T?(Q 2 ) = C NS (Q 2 ) (f| + ||) + \c s {Q 2 ) (E + 2Ac) 



(2) 



where Ac is 

s M Ac<y ) = (p, s I (c7 m 7 5 c) m2 I p, s) (3) 

at a relevant mass scale /i 2 . Here the coefficients Cns{Q 2 )i Cs(Q 2 ) account for pertur- 
bative corrections of the form 1 — a s /n + . . . 0. The measureable singlet axial coupling 
constant can be defined as 

gf = AS + 2Ac (4) 

where the first term is the contribution of all light quarks (u,d,s) together, defined as 
AS = An + Ad + As. 

Both charm and light quark contributions can be calculated from the divergence of 
the axial current: for each flavor one has (in Euclidean space) 

- dniqEltilhQE) = 2m q q E jsqE + T7r^K G ^^)E (5) 

l07T z 

However, the role of these two terms is different for each flavor. For light quarks one 
can take the chiral limit m — ► and ignore the first term in the r.h.s. The OPE 



1 See e.g. || for a recent review and nomenclature. 
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expansion of the first term in 1/m was done by Halperin and Zhitnitsky || and it 
shows that the anomaly term actually cancels, while the leading term in the heavy 
quark limit is 

- d,(c El ^c E ) = +l —^f^G%G b ua G c au + 0(G 4 /m 4 c ) (6) 

Note that it is indeed suppressed by l/m 2 , as it should. 

3. Using lattice gauge configurations (or a model which defines gauge fields in some 
way, as we do below), one can correlate these two gluonic operators with the nucleon 
and thus evaluate the matrix elements in question, AE and Ac. Information about 
them is contained in the three point correlation functions 

n± {u} (-x/2,x/2,y) = T% }± (0 | T{ Va (-x/2)j^(y)rj^(x/2)} | 0) (7) 

where \i is the component of the axial current j^(y) = qi^l^q and {v} a multiple n- 
index for the projection operator = 7^1 . . . r ) vn P±- l P± = Pr±Pl- The r) a (x), rf(x) 
are the so called Ioffe currents, local operators with the quantum numbers of the 
nucleon or in a simplified consideration the quarks. 

Before we go for the evaluation of the three-point functions, let us qualitatively 
explain the signs and magnitude of the gluonic operators for an instanton (anti- 
instanton) configuration. Introducing the topological charge Qi=l (-1) and substi- 
tuting the known analytic solution for the instantons one finds 

J-r g -o 1 12p4 m 



and 



1 ~ Qfi 1 (i 6 

£abcr~ia r~ib r~ic s-\ f_ (Q\ 



2 Since the latter do not contain charm quarks, the contribution of charm quarks is only via 
the disconnected, OZI violating, diagrams. 

4 



where p is the size and the zi the position. The first important observation now is that 
each instanton/anti-instanton contributes to AS and Ac with the opposite sign. 

The second is that, by dimension, charm component relative to light quark one is 
~< p~ 2 > /m 2 , where angular brackets mean averaging over instanton size distribu- 
tion. Obviously this average favor smaller sizes than the usual r.m.s. one. 

Of course, due to the factor of the topological charge Qi = ±1 the average of 
both of those pseudo scalar operators in the vacuum vanishes. Non-zero results can 
only be obtained if it is correlated with another pseudo scalar observable, such as aq, a 
product of spin vector and the momentum transfer. Another useful way to explain it is 
to note that a spin vector acts on the vacuum like a dipole operator, thereby polarizing 
(topological) charges and creating non-equal densities of instantons and anti-instantons 
along its vector. 

4. Quantitative calculations are based on the three-point correlation functions indi- 
cated above. Instead of using the axial current itself, we are considering its divergence 
(|5|) in eq. (|7|) . Then for x in Euclidean time direction, the spin asymmetry induced by 
the pseudoscalar gluonic operators in the direction n of y can be traced out with T™ - . 
The correlator is then ~ (y ■ s) which translates into ~ (q ■ s) in momentum space. 
Since we plot the correlator in coordinate space as a function of y we are not facing 
the problem of a vanishing forward matrix element J7| of these operators. 

For comparison, we use two instanton ensembles, the random instanton liquid 
(RILM) and the so called interacting instanton liquid model (IILM) (for review see 
||). The former is a simple model based on phenomenological parameters of the in- 
stantons like the average size, while the latter is a sophisticated model including the 
't Hooft effective interaction to all orders. The ensembles have 256 instantons in a 
rectangular box of (5.7fm) 2 x (2.8fm) 2 . The instanton density is the same n = 1 fm~ 
for all configurations. 
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Correlations functions with the nucleon propagator decay very fast with x due to 
large nucleon mass, which in practice strongly restricts the range of x,y at which the 
measurements can be made. Virtual charm pairs are rather well localized anyway, 
and so one may probably use a single "constituent" quarkQ instead of three in the 
nucleon as a reasonable first approximation. That is what we did below. So we do 
not actually discuss nucleons, but some "generic" hadrons, made of loosely bound 
constituent quarks. 

Our results for both correlators are shown in Fig. 1 and Fig. 2. The geometry of 
the x and y directions and their ratio we have adopted from previous three-point corre- 
lation function studies [ Hj : their particular choice is not important for the conclusions. 



The dashed lines are guiding analytic calculations based on the single instanton ap- 
proximation, and serve only as a benchmark for the numerical data at smaller distances 
(x, y <C 0.5fm). 

The main finding is that both operators show very similar shape of the correlation 
functions. This is further illustrated in Fig.3, where we plot their ratio. The instanton 
ensemble used however matters: for IILM the charm-related operator shows a much 
weaker signal. This is related to different instanton size distributions: in RILM all 
of them have p = l/3fm while in IILM the distribution is peaked at about 0.4fm. 
Another important reason for the reduction effect in IILM are positive (screening-type) 
correlations between instantons and anti-instantons induced by quark exchanges. We 
use the similarity in shape of both correlators as a sign that the form factors for both 
S(g 2 ) and Ac(q 2 ) are the same, and the ratio of the correlation functions shown in 
Fig.3 can be immediately translated to the ratio ^p. 



3 More complete analysis of structure functions of a constituent quark resulting from those 
measurements will be presented elsewhere M. 
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Finally our results are 



Ac 

— ~ -0.20 ± 0.04 (RILM), (10) 

Zj 

and 

Ac 

— ~ -0.08 ±0.006 (IILM). (11) 

Zj 

Several things have to be remarked now. Using an approximate scheme for the 
non-zero mode propagator involved in the correlation function would give the analytic 
expression ^ = — 57^7 {mp) 2 — — 0-2 |§ for the standard values of the RILM. This 
scheme is also shown in Fig. 3 and actually very close to the full RILM simulation. 

The positive sign of AS here agrees with the experimental data and the negative 
sign for Ac is somehow in agreement with lattice calculations [11] of disconnected 



matrix element of the axial current itself, which give a negative sign for the three light 
quark flavors. Both AS and Ac signs are opposite however to the scenario proposed 
by Halperin and Zhitnitsky ||, which is mainly based on an older low energy theorem 
by Kiihn and Zakharov |12| . 



These numbers ((U^O) are significantly smaller than the crude estimates for gy 
based on Goldberger-Treiman-type relations for rf exchange used in Nevertheless 
the size of Ac is in fact surprisingly large. It seems to be only 3-6 times smaller 
than that of the strange quark sea. Hopefully the next generation of polarized DIS 



experiments [13,14] with charm-jet tagging will be able to observe it. 
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Fig. 1: The quark three-point correlation function for the charm contribution ~ 
f abc G a (y)G b (y)G c (y) for the random configuration (RILM). The approximate (app) for- 
mula is also shown. 
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Fig. 2: The quark three-point correlation function for the charm contribution ~ 
f abc G a (y)G b (y)G c (y) for the interacting configuration (IILM). 
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Fig. 3: The ratio of the three-point function for ^?GG and —^f abc G a {y)G b {y)G c {y) 
for the RILM and IILM configuration. 
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